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1.  INTRODUCTION  AND  SUMMARY 

This  report  discusses  the  effects  that  large  amounts  of 
clay  and  water  in  the  surrounding  rock  may  have  on  the  contain- 
ment of  hot  cavity  gases  due  to  an  underground  nuclear  event. 

The  containment  failure  of  the  Baneberry  event,  detonated  in 
Yucca  Flat  in  December  of  1970,  seems  to  be  an  ideal  vehicle 
for  this  study.  The  Baneberry  site  has  been  characterized  by 
the  USGS^  as  having  unusually  high  concentrations  of  mont- 
morillonite  clay  (65  percent)  around  the  working  point.  The 
water  content  of  the  working  point  material  (an  altered  tuff) 
is  related  to  the  clay  content  and  has  been  estimated  by  both 
the  USGS  ^ and  by  L.  Ramspott^  to  be  as  much  as  25  percent 
by  weight.  Material  properties  data  for  Rainier  Mesa  tuff, 
compiled  by  S.  Butters  of  TerraTek,  Inc. ^ indicate  that  high 
water  content  in  laboratory  samples  correlates  with  low  shear 
strength.  While  there  is  little  hard  data  on  the  effect  of 
clay  on  material  strength,  it  seems  clear  that  the  presence 
of  clay  can  only  weaken  the  host  material. 

A scenario  for  the  containment  failure  of  Baneberry  has 
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been  presented  in  a report  by  N.  Rimer  which  discusses  one- 
dimensional ground  motion  calculations  using  the  SKIPPER  code 
of  the  Baneberry  event  and  of  typical  Rainier  Mesa  and  Yucca 
Flat  events  which  were  contained.  Material  properties  for 
Baneberry  were  based  on  site  investigations  and  analyses  by 
USGS^  and  L.  Ramspott^  which  led  to  the  assumption  of  three 
spherical  layers  about  the  Baneberry  working  point.  The  working 
point  layer  of  altered  tuff  was  assumed  to  have  negligible  shear 
strength  due  to  the  very  high  clay  and  water  content.  The  second 
and  third  layers  for  the  calculations  were  both  alluvium  (the 
third  layer  was  of  infinite  extent) . Since  the  clay  content  of 
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the  second  layer  had  been  estimated  at  between  20  and  50  per- 
cent, it  was  assumed  to  have  half  the  shear  strength  of  the 
third  layer. 

The  major  difference  between  the  Baneberry  calculation 
and  the  Rainier  Mesa  and  Yucca  Flat  results  appeared  to  be  in 
the  magnitude  of  the  compressive  residual  stress  fields  due  to 
nonuniform  plastic  loading  and  subsequent  unloading  of  the  rock 
around  the  explosion-produced  cavities.  Residual  stress  fields 
are  seen  after  the  cavity  rebounds  in  all  ground  motion  calcu- 
lations at  S3.  Parameter  studies  by  N.  Rimer ^ have  indicated 
that  the  peak  transverse  residual  stresses  depend  very  strongly 
on  the  shear  strength  of  the  rock.  Therefore,  it  was  not  sur- 
prising that  the  Baneberry  calculation,  modeled  with  relatively 
low  strength,  gave  significantly  lower  residual  stresses. 

Based  on  these  results,  the  following  scenario  was  pre- 
sented as  possibly  leading  to  the  venting  of  cavity  gases  to 
the  surface  at  late  times.  The  low  strength  due  to  high  clay 
and  water  content  led  to  low  residual  stress  fields.  In  the 
absence  of  a significant  residual  stress  field,  the  hot  cavity 
gases  could  cause  tensile  hydrofracturing.  The  fracturing, 
if  sufficiently  great,  could  lead  to  venting.  The  low  material 
strength  was  considered  to  be  the  primary  reason  for  the  con- 
tainment failure  of  Baneberry. 

Although  the  original  report  of  this  work  was  written  in 
the  summer  of  1975,  it  received  only  limited  distribution  to 
interested  ERDA  representatives  at  the  time.  Because  of  the 
general  sensitivity  of  the  subject  matter,  we  were  asked  to 
withhold  further  distribution  until  additional  calculations  were 
completed  which  would  determine  the  sensitivity  of  the  results 
to  the  modeling  assumptions.  A set  of  suggested  multi-layer 
calculations  was  completed  and  presented  at  an  Earth  Motion 
Calculators  meeting  held  in  La  Jolla  on  May  4,  1976.  Since  the 
multi-dimensional  calculations  requested  in  Ref.  4 have  recently 


been  done  and  presented  by  LLL  at  the  66th  CEP  meeting,  the 
report  was  released  in  January  of  1977. 

In  the  intervening  period,  considerable  work  has  been 
done  to  clarify  questions  raised  by  that  report.  The  present 
work  includes  a discussion  of  constitutive  modeling  and 
material  properties  data  improvements  (see  Section  2)  used  in 
later  calculations  of  both  Baneberry  and  the  comparison  events 
in  Rainier  Mesa  and  Yucca  Flat.  It  was  felt  that  the  com- 
parisons in  Ref.  4 might  not  be  completely  valid  since  there 
were  some  modeling  differences  in  the  calculations  being  com- 
pared (the  Rainier  Mesa  and  Yucca  Flat  calculations , as  des- 
cribed in  Ref.  4,  were  not  made  specifically  as  part  of  this 
program.)  Careful  consideration  has  shown  that  those  comparisons 
are  qualitatively  valid. 

Section  3 describes  the  results  of  calculations  presen- 
ted at  the  Earth  Motion  Calculators  meeting  of  May  4,  1976. 

Based  on  suggestions  from  ERDA  representatives,  the  sensitivity 
of  the  results  to  some  of  the  modeling  assumptions  used  in  the 
Baneberry  calculation  was  investigated.  It  was  found  that  the 
basic  results  of  Ref.  4 were  not  sensitive  to  any  parameter 
other  than  material  strength. 

In  Section  4,  calculations  are  presented  for  a weak  and 
a strong  Rainier  Mesa  saturated  tuff,  a weak  and  a strong 
Yucca  Flat  dry  tuff,  as  well  as  a weak,  best-guess  Baneberry 
calculation  and  an  upper-limit,  relatively-high-strength  Bane- 
berry calculation.  These  calculations  show  Baneberry  to  have 
a much  lower  residual  stress  field  than  the  other  NTS  events 
due  to  the  lower  strength  of  the  surrounding  alluvium. 

Section  5 discusses  "figures  of  merit,"  measures  by 
which  one  can  tell  whether  a particular  event  is  a greater 
containment  risk  than  other  previously-detonated  successful 
events.  The  figure  of  merit  chosen  was  the  ratio  of  the  cavity 


5 


pressure  to  the  peak  transverse  residual  stress.  For  this 
comparison,  a sophisticated  determination  of  the  cavity 
pressure  was  used  as  described  below.  Using  the  resulting 
figure  of  merit,  Baneberry  was  clearly  the  greatest  contain- 
ment risk  of  the  events  examined. 


2.  CONSTITUTIVE  MODELING  AND  MATERIAL  PROPERTIES 


In  Section  2.1,  the  constitutive  models  are  discussed 
briefly  with  emphasis  on  the  modeling  improvements  between 
Ref.  4 and  the  calculations  presented  here.  Section  2.2  gives 
the  material  properties  used  for  the  calculations  of  Ref.  4 
and  of  Section  3 of  this  report.  The  modeling  and  material 
properties  used  for  the  results  presented  in  Section  4 are 
discussed  in  Section  2.3. 


2.1  CONSTITUTIVE  MODELING 

Descriptions  are  given  here  of  the  constitutive  models 
used  in  the  SKIPPER  calculations.  A more  complete  discussion 
may  be  found  in  Cherry,  et.  al46^ 


Equations  of  State 

The  Baneberry  calculations  described  in  Ref.  4 used  an 
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equation  of  state  of  the  form  developed  by  Tillotson  and 
fit  at  S3  to  available  shock  lata  for  a number  of  geological 
materials.  For  compressed  states  (p  > p ) and  for  cold  expan- 
ded states  (p  < pQ  and  e < eg) , the  pressure  is  given  by 


P 

s 


a + 


b 


— + 1 

V1* 


ep  + Ay  + By2 


For  expanded  states  (p  < pQ)  where  e > e'  , the  pressure  is 
given  by 


P 

\r 


aep 


4 


bep 


+ Ay  e 


The  phase  transition  from  liquid  to  vapor  (p  < pQ  and 

e < e < e')  is  approximated  as 
s s 

P = T(e  - e ) P + (e'  - e)  P 1 . 

eg  ~ eg  L s v s s J 

Here 

p = mass  density 
n = p/pq 
y = n + l 

eg  = specific  internal  energy  as  the  material  is 
brought  to  vaporization  temperature 

e'  = additional  specific  internal  energy  required 
to  change  the  material  from  a liquid  to  a 
vapor  state 

All  the . above  formulas  assume  that  no  air-filled  voids  are 
present  in  the  rock.  The  constants,  a,  b,  and  eQ  govern  the 
energy  dependence  of  the  material,  allowing  it  to  behave  as 
an  ideal,  y-law  gas  at  high  energies  (a  ^ = y-1) , and  reducing 
the  effective  y at  lower  energies.  The  density  dependence  is 
governed  by  a quadratic  polynomial  with  constants  A and  B 
which  may  be  fit  to  shock  data.  The  Tillotson  equation  of 
state  does  not  accurately  model  the  expansion  from  rock  states 
above  approximately  10  kbars  where  the  expansion  of  steam  is 
important. 

The  Tabular  Array  of  Mixtures  Equation  of  State  (TAMEOS) 
described  in  Cherry^6 ^ has  been  used  for  all  tuff  calculations 
discussed  here.  TAMEOS  mixes  the  proper  percentage  of  water 
with  grain  density  rock,  assuming  pressure  equilibrium  between 
the  rock  and  the  water,  to  give  the  pressure  response  for  the 
fully  crushed  tuff  mixture.  This  gives  a far  better  description 
of  the  subsequent  expansion  of  the  shocked  water-rock  mixture. 
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Porous  Crushup 


The  equations  of  state  discussed  above  describe  the 
pressure  response  for  fully  saturated  rocks.  Air-filled 
porosity  <f>0  is  included  through  .the  S3  porous  crushup  model 
(P-a  model)  in  which  the  pressure  of  the  porous  material 
is  described  by 


where 

V = specific  volume  of  porous  material 

e = specific  internal  energy  of  porous  material 
P = pressure  obtained  from  equation  of  state 
a = distension  ratio  defined  by  V/V  >_  1 

V = specific  value  of  the  material  with  zero  air- 

filled  voids 

The  distension  ratio  is  required  to  decrease  from  an 
initial  value 


at  zero  pressure  down  to  1.0  as  the  pressure  increases  to  P , 
the  crush  pressure  or  pressure  limit  at  which  all  air-filled 
porosity  is  assumed  irreversibly  removed.  Any  porosity  lost 
during  compression  below  the  elastic  pressure  Pg  is  recover- 
able upon  expansion.  Any  porous  crushup  occuring  between  Pg 

and  P is  irreversible, 
c 

Elastic  Properties  and  Material  Strength 

The  ambient  compressional  wave  velocity  CQ  is  related 

to  the  bulk  modulus  K and  the  shear  modulus  G by  the 

o 

relationship 
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2 4 

p C = K + t G 

'o  o o 3 

Here,  the  shear  modulus  is  assumed  constant.  Since  the  bulk 
modulus  K varies  with  pressure,  Poisson's  ratio  a is  not 
being  held  constant.  One  calculation  was  made  with  constant 
a throughout  the  crushup.  This  requires  a constant  ratio  of 
K/G,  i.e.,  a varying  G. 

The  deviatoric  components  of  the  stress  tensor  are 

related  to  the  deviatoric  strain  rates  by  Hooke's  law 

Sij  = 2 G e.. 

The  material  strength  model  requires  that  the  principal  stress 
be  within  the  von  Mises  yield  surface.  For  spherical  geometry, 
this  reduces  to  the  condition  that  the  magnitude  of  the  radial 
deviatoric  stress  not  exceed  2/3  the  yield  strength  (maximum 
stress  difference) . Behavior  of  a failed  element  is  governed 
by  the  non-associated  flow  rule. 

For  the  calculations  of  Ref.  4,  the  yield  strength  was 
given  by 

Y = [y  + Y 5—  (2  - Y|  (l  - — ) P < P , e < e 

L 0 m Pm  ' V-l  ' em  ’ m m 

v ■ (*. + ht)  P-Ve<  % 


Y = 0 


e > e 


where  Y , Y , P and  e are  constants  for  a given  material, 
o m m m 

The  more  recent  calculations  use  an  improved  form  for 
the  yield  strength  which  has  been  shown  by  Cherry  and  Peterson^1 
to  be  empirically  superior  to  a simple  pressure  dependence. 
Here,  the  pressure  P in  the  above  equations  is  replaced  by 
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P which  is  a function  of  P and  the  deviatoric  stresses 


given  by 


P = P 


-M 


S11S22S33 


A description  of  how  shear  failure  is  treated  in  this  model  is 
given  in  Cherry,  et.  al_.  ^ 

Tension  failure  is  allowed  to  occur  in  an  element  if  a 

principal  stress  becomes  tensile.  The  tension  failure  model 
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proposed  by  Maenchen  and  Sack  and  described  in  detail  by 
Cherry^  is  then  applied.  An  inelastic  strain  which  is  just 
sufficient  to  zero  the  tensile  stress  is  introduced  normal  to 
the  crack.  This  strain  increases  or  decreases  as  the  crack 
opens  or  closes. 


Cavity  Equation  of  State 


The  cavity  source  region  was  chosen  to  be  initially 
large  enough  to  vaporize  70  metric  tons  of  rock  for  each  kilo- 
ton  of  device  yield.  (The  rock  is  assumed  to  be  at  undisturbed 
density.)  The  cavity  pressure  for  all  calculations  is  given  by 

P = (y-1)  e/V 

For  the  Rainier  Mesa  and  Yucca  Flat  calculations  des- 
cribed in  Ref.  4,  y was  set  equal  to  1.5.  The  Baneberry  cal- 
culation of  Ref.  4 assumed  that  y was  a function  of  specific 
volume  given  by 


y = 1.03  + 0.9//V 


This  expression  for  y was  a preliminary  result  obtained  for 
tuff  modeled  using  the  chemical  equilibrium  CHEST  code  oi 


D.  Laird 


and  has  been  modified  for  later  calculations.  It 


will  be  shown  later  that  the  peak  residual  stresses  are  not 
sensitive  to  the  choice  of  cavity  y.  However,  depending  upon 
the  cavity  model  used,  peak  velocities  at  a given  location  may 
differ  by  20  percent. 

More  recent  calculations  used  the  following  expression 

for  y • 


y'  = 1.085 

y'  = 1.2085  - 0.065  V 
y'  = 1.9375  - 0.875  V 
y = y'  - 0.032  + 3.2  x 10_13  e 
y = y' 


V > 1.9 

1.9  > V >_  0.9 

V < 0.9 

e > 1011  ergs/gm 
e < 1011 


It  should  be  emphasized  that  the  above  expressions  for 
Y are  only  very  approximate  fits  to  the  complicated  tabular 
equation  of  state  of  Laird.  The  fit  given  above  was  derived 
from  data  from  a table  developed  for  one  particular  water 
content  (24  percent  by  weight,  including  approximately  7 per- 
cent bound  water) . Since  the  creation  of  a CHEST  table  is  a 
complicated  procedure,  the  same  expression  for  y was  used  in 
calculations  for  Baneberry,  Rainier  Mesa  and  Yucca  Flat  tuffs, 
i.e.,  over  a wide  range  of  water  contents. 


2.2  MATERIAL  PROPERTIES  FOR  THE  PARAMETER  STUDY 

Table  1 gives  the  material  properties  data  used  in  the 
original  calculations  of  Ref.  4.  For  Baneberry,  Material  1 
describes  a saturated,  high  clay  content  altered  tuff  extend- 
ing radially  15.24  meters  (50  feet)  from  the  working  point. 
Since  the  clay  content  of  this  layer  was  greater  than  50  per- 
cent, the  material  strength  was  considered  negligible.  The 
Tillotson  equation  of  state,  using  constants  approximating  the 
data  of  Stephens,  et . al^.  , was  used  to  model  the  pressure 

response  of  the  material. 


12 


rn  t> 
<u 
to  o 
O 3 
O TJ 
3 O 
P 
Q- 

- a) 

P 

P 

<H  <D  <P 
X)  <P 
a)  aj  3 

i— I C 4J 

X»  <0 
3 CQ  "O 
Eh  <D 
P -P 
O <0 
>P  P 
3 
(0  4-> 

aj  « 

•H  W 

x> 

p to 

Q)  V) 

a a) 
o S3 
p 

a P 

a> 


P 

dl  <0 

rH 

■*r 

CM 

rH 

03 

in 

in 

in 

VO 

in 

•H  W 

a* 

co 

r^» 

CM 

O 

rH 

o 

^3* 

o 

in 

o 

o 

rH 

c <u 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 1 1 1 1 

i i i i 

Tj* 

« 

rH 

CM 

35 

57 

o 

O 

o 

o 

o 

rH 

o 

CM  I 1 1 1 

1 1 1 1 o 

<0  -P 

13 

00 

in 

r-H 

62 

12 

18 

36 

Os 

O 

o 

ro 

rH 

O d) 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 1 1 

i i i 

i i 

U f-l 

3 (P 

>H 

rH 

CM 

25 

55 

o 

o 

o 

O 

o 

o 

CM  | | 

i i i 

1 1 o 

o 

O O 

o 

r* 

VO 

in 

VO 

rH  rH 

o 

o 

o 

o 

rH 

r- 

in 

ro 

rj* 

vo 

rH 

in 

CM  O O 

o 

o 

o 

CM 

o 

CM 

ro 

r- 

O 

a» 

ro 

o 

o 

rH 

ro 

rH 

rH 

r- 

o 

00 

VO 

o 

00 

in  in 

o 

r- 

VO 

in 

O 

o o 

o 

r- 

o 

O 

o 

rH 

in 

ro 

H* 

VO 

CM 

o 

o o o 

o 

o 

o 

CM 

o 

CM 

ro 

O 

o> 

ro 

O 

o 

ro 

in 

rH 

rH 

r- 

o 

00 

o 

o 

VO 

in 

VO 

o 

rH 

in 

vo 

• 

i 

1 

• 

• 

• 

• 

• 

• 

rH 

o o o 

O 

o 

i 

1 

o 

CM 

ro 

O 

o 

CM 

o 

o 

** 

rH 

rH 

in 

ro 

§»  i 

g 

O’ 

I 

\ 

\ 

\ 

\ 

tP 

tP 

CP 

CP 

* 

p 

p 

p 

p 

o 

<u 

di 

d) 

d> 

0) 

. s 

w 

^-S 

p 

p 

p 

p 

p 

P 

p 

o 

o 

O 

o 

"s 

P 

ctf 

CT3 

<a 

<TJ 

<d 

rH 

rH 

rH 

rH 

p 

^ — . 

g 

<TJ 

XI 

XI 

XI 

XI 

XI 

XI 

XI 

O 

O 

O 

o 

(TJ 

p 

XI 

-SC 

-SC 

x 

-V 

rH 

rH 

rH 

rH 

X 

<TJ 

X 

s— 

v_^ 

■H— ^ 

XI 

.* 

o 

0 

UI 

o 

0 

g 

d> 

O 

c 

> g 

0)  ' W o 

— 

o 

o 

k 

Y 

>H 

a 

a 

& 

0< 

at 

d> 

<u 

di 

< 

CO 

<0 

XI 

Materials  2 and  3 are  assumed  to  be  alluvium  modeled 
using  a Tillotson  equation  of  state.  Material  2,  with  greater 
than  20  percent  montmorillonite , is  assumed  to  have  half  the 
shear  strength  of  Material  3,  which  is  a weak  alluvium. 

Material  2 was  assumed  to  be  36.58  meters  thick  (120  feet) 
and  to  be  surrounded  by  Material  3.  Free  surface  effects  and 
a radially  varying  gravitational  field  were  neglected  in  these 
calculations. 

[1  21 

The  measured  data  ' upon  which  the  modeling  was  based 
was  often  conflicting  since  the  Baneberry  geology  is  unusually 
complicated  and  since  the  data  were  accumulated  from  holes 
drilled  both  pre  and  postshot  and  in  different  locations  rela- 
tive to  the  WP.  The  log  data  indicated  that  Material  2 had  a 
compressional  wave  velocity  20  to  30  percent  greater  than 
Materials  1 or  3.  This  was  included  in  the  modeling.  Also, 
due  to  the  higher  clay  content,  Material  2 had  only  half  the 
air-filled  void  content  of  Material  3.  Also  included  in  Table  1 
is  material  properties  data  used  for  calculations  for  Yucca 
Flat  dry  tuff  and  for  Rainier  Mesa  saturated  tuff.  These  cal- 
culations were  not  made  specifically  for  the  Baneberry  study 
but  were  used  for  comparison  with  the  Baneberry  results.  As  a 
result,  differences  exist  in  the  constitutive  modeling  for  the 
three  calculations.  For  example,  both  Yucca  Flat  and  Rainier 
Mesa  calculations  use  a TAMEOS  equation  of  state  rather  than  a 
Tillotson  equation  of  state. 

A parameter  study  was  made  to  examine  the  sensitivity  of 
the  results  of  Ref.  4 to  the  modeling  assumptions  used  in  the 
Baneberry  calculation.  The  starting  point  for  this  study  is 
the  material  properties  data  of  Table  1.  The  calculations  of 
Ref.  4 all  used  the  pressure  dependent  failure  surface  described 
in  Section  2.1.  However,  the  parameter  study  used  the  empiri- 
cally superior  P dependent  failure  surface.  The  parameter  study 
also  used  the  newer  functional  dependence  of  y on  specific  volume 
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and  energy  which  was  described  in  Section  2.1  for  the  cavity 
equation  of  state. 

2.3  MATERIAL  PROPERTIES  FOR  CALCULATIONS  OF  SECTION  4 

Using  the  constitutive  modeling  described  in  Section  2.1, 
a series  of  calculations  was  made  for  Rainier  Mesa  (Dido  Queen 
and  Mighty  Epic  events) , for  Yucca  Flat  dry  tuff  (a  calcula- 
tion using  average  properties)  and  for  Baneberry.  Table  2 
gives  the  properties  used  for  these  calculations.  All  calcu- 
lations use  the  TAMEOS  equation  of  state  for  the  working  point 
material,  a cavity  equation  of  state  where  y is  a function 
of  specific  volume  and  energy,  and  the  P dependence  for  the 
failure  surface.  This  section  discusses  the  properties  chosen 
for  each  event. 

The  material  properties  for  the  Dido  Queen  event,  den- 
sity, grain  density,  water  content  and  compressional  wave 
velocity  were  obtained  from  the  CEP  document  for  that  event. 

The  failure  envelope  was  modeled  using  Figure  24  of  Butters, 

f 12 1 

et . al . while  the  crush  curve  was  obtained  from  Figure  19 

of  the  same  report.  Dido  Queen  was  chosen  for  this  study 
because  it  was  located  in  tunnel  e of  the  Rainier  Mesa  complex, 
a tunnel  that  triaxial  loading  tests  indicate  contains  rela- 
tively weaker  tuffs. 

The  Mighty  Epic  event,  located  in  tunnel  n,  was  chosen 
as  an  example  of  a stronger  Rainier  Mesa  tuff.  The  choice  of 
material  properties  for  the  Mighty  Epic  event  is  discussed  in 
detail  in  Rimer,  et.  al. 

The  Yucca  Flat  material  properties  were  based  on  average 

[14] 

properties  for  a Yucca  Flat  dry  tuff  summarized  by  Ramspott. 
Unfortunately,  there  is  no  data  on  the  strength  of  this  mate- 
rial. TerraTek  laboratory  data  for  Rainier  Mesa  tunnel  bed 
tuffs  compiled  by  S.  Butters ^ shows  a relationship  between 
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shear  strength  and  water  content.  It  should  be  emphasized 
that  the  data  of  Ref.  3 is  stress  difference  in  uniaxial  strain 
at  4 kbars  confining  pressure,  not  the  yield  strength.  However, 
it  should  give  a relative  indicator  of  the  strength  of  the 
materials . 

F.  App^1^  has  drawn  several  straight  lines  through  a 
plot  of  strength  (stress  difference  of  4 kbars)  v£  water  con- 
tent (by  weight)  for  these  data  which  indicate  that  the  shear 
strength  is  inversely  proportional  to  the  water  content. 

These  lines  shown  in  Figure  1 imply  that  for  Baneberry  water 
contents,  the  strength  is  quite  small.  For  the  lower  water 
content  of  Yucca  Flat  tuff  above  the  water  table,  the  plot 
gave  a strength  value  of  0.63  kb.  The  yield  strength  of 
0.63  for  Yucca  Flat  tuff  water  content  was  apportioned  between 
Yq  and  Y^  in  the  ratio  suggested  by  the  Dido  Queen  data.  This 
leaves  the  magnitude  of  yet  to  be  determined.  For  this 

reason,  two  calculations  were  made  for  Yucca  Flat,  one  for  a 

P of  3 . 0 kbar,  as  in  Dido  Queen,  and  another  for  a P of  1.0 
m m 

kbar,  as  in  Mighty  Epic.  The  parameter  Pm  is  as  important  to 
the  residual  stress  field  as  the  magnitude  of  the  failure 
strength  itself.  A low  value  of  Pm  indicates  a steeper  slope 
of  the  failure  surface  and  therefore  a higher  strength  at  low 
pressures. 

A Yucca  Flat  calculation  made  with  the  data  of  Table  2 
using  the  TAMEOS  equation  of  state  was  compared  to  one  made 
using  the  Tillotson  equation  of  state.  Tillotson  constants 
A and  B for  this  calculation  were  chosen  to  match  the  zero 
pressure  bulk  modulus  used  in  the  TAMEOS  calculation  and  to 
match  its  zero  energy  curve  at  the  crush  pressure  P . The 
comparisons  are  given  in  Table  3. 
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Stress  Difference  (bars) 


Table  3 


Comparison  of  results  using  TAMEOS  and  Tillotson 
equations  of  state  for  Yucca  Plat  dry  tuff 


Cavity 

radius 

(meters) 

Cavity 

pressure 

(bars) 

Radial 

stress 

(bars) 

Transverse 

stress 

(bars) 

TAMEOS 

31.70 

136.7 

174 

304 

Tillotson 

29.18 

183.9 

140 

249 

The  calculations  showed  about  a 10  percent  difference 
in  cavity  radius  and  close  to  20  percent  differences  in 
radial  and  transverse  peak  residual  stresses.  Cavity  pres- 
sures were  even  further  apart.  On  the  basis  of  these  results, 
it  was  decided  to  apply  TAMEOS  to  the  Baneberry  working  point 
material.  Input  to  TAMEOS  is  the  water  content  (by  weight) , 
the  grain  density  of  the  tuff,  and  constants  defining  the 
equation  of  state  of  the  grain  density  tuff.  It  was  assumed 
that  the  equation  of  state  constants  applicable  to  Rainier  Mesa 
tuff  (having  much  lower  grain  densities)  may  be  used  for  the 
Baneberry  altered  tuff.  The  important  features  are  the  high 
water  content  and  the  high  grain  density  due  to  clay  content. 

Since  the  peak  residual  stresses  had  been  shown  to  be 

primarily  dependent  on  the  failure  strength  and  not  to  be  very 

dependent  upon  the  exact  position  of  the  layer  interface  or  on 

the  compressional  wave  velocity  (the  results  are  presented  in 

Section  3) , it  was  decided  to  model  the  Baneberry  site  using 

only  two  layers,  a working  point  layer  of  over  50  percent 

montmoril Ionite  extending  radially  35.66  meters  (117  feet)  from 

the  WP , with  an  alluvium  layer  of  lower,  but  still  significant 

clay  content  around  it.  The  location  of  the  layer  was  based 

r 2] 

on  data  from  drill  hole  Ue8i  discussed  by  Ramspott.  The 

Tillotson  equation  of  state  was  used  for  the  second  layer  which 
is  not  expected  to  see  stresses  greater  than  5-10  kbars.  Both 
layers  were  modeled  using  data  from  hole  Ue8i. 


19 


I -i . 


For  layer  2,  data  on  the  failure  envelope  has  been 
found  for  Diagonal  Line  alluvium,  considered  to  be  a rela- 
tively strong  alluvium  and  for  Merlin  alluvium,  consiciered 

to  be  a more  representative  alluvium.  The  strength  given  for 
layer  2 for  the  calculation  labeled  Baneberry  1 in  Table  2 
was  obtained  using  Diagonal  Line  alluvium  data  and  so  repre- 
sents an  upper  limit  to  the  strength  of  layer  2.  The  calcu- 
lation labeled  Baneberry  2 is  identical  with  Baneberry  1 
except  that  the  strength  of  layer  2 is  lower  as  shown  since 
it  was  modeled  using  Merlin  data. 

The  failure  strength  of  the  Baneberry  working  point 
layer  was  modeled  using  data  from  Stephens,  et.  al. 

Figure  7 of  Ref.  17  plots  shear  strength  v£  water  contents  for 
alluvium.  That  plot  indicates  that  for  water  contents  of  20 
percent  or  more,  alluvium  has  zero  shear  strength.  This  implies 
that  the  working  point  material  may  have  negligible  strength. 
Figure  1 indicates  negligible  strength  for  water  contents  of 
30  percent. 


3.  SENSITIVITY  OF  BANEBERRY  RESULTS  TO  MODELING  ASSUMPTIONS 


The  results  of  a parameter  study  presented  at  the 
Earth  Motion  Calculators  meeting  held  in  La  Jolla  on  May  4, 
1976  are  discussed  in  this  section.  The  purpose  of  this 
study  was  to  answer  questions  raised  by  ERDA  representatives 
as  to  the  sensitivity  of  the  results  of  Reference  4 to  the 
modeling  assumptions  used  for  the  Baneberry  calculation.  The 
material  properties  used  for  the  3 material  Baneberry  base 
calculation  are  those  used  for  Reference  4 and  are  given  in 
Table  1 of  Section  2.2.  A second  look  at  the  modeling 
assumptions  indicated  that  the  thickness  of  layer  2 (mate- 
rial 2)  should  be  51.8  meters  (170  feet)  rather  than  36.58 
meters  (120  feet)  as  in  Reference  4.  Modeling  improvements 
were  incorporated  into  this  study  as  described  in  Section  2. 

\_  The  calculations  of  this  study  are  the  following: 

1.  The  base  calculation:  This  calculation  uses  the  data 

of  Table  1 but  assumes  that  layer  2 is  170  feet  thick.  An 
important  feature  is  that  the  cavity  gas  equation  of  state 
has  a constant  y of  1.5. 

The  yield  strengths  are 


Layer 

1 

(clay) 

o 

li 

Layer 

2 

(half  strength) 

Y =7.0 
o 

bars. 

Y =50  bars, 
m 

P = 700 
m 

bars 

Layer 

3 

(full  strength) 

Y =7.0 
o 

bars. 

Y = 100  bars, 
m 

Pm  = 700 
m 

bars 

la.  This  calculation  is  identical  with  the  base  calcula- 
tion except  that  layer  2 is  120  feet  thick.  It  essentially 
duplicates  the  Baneberry  calculation  of  Reference  4 with 
better  modeling  for  the  failure  surface  and  with  a constant 
y for  the  cavity  gas. 


21 


2.  High  strength,  constant  y:  This  calculation  is 

identical  with  calculation  (1)  except  for  material  strengths. 
The  new  strengths  are 


Layer 

1 

Y = 
o 

30 

bars, 

Ym  = 10 

bars. 

P = 700 
m 

bars 

Layer 

2 

It 

0 

>< 

30 

bars. 

Ym  " 50 

bars. 

P = 700 
m 

bars 

Layer 

3 

n 

o 

>1 

30 

bars. 

Ym  " 100 

bars 

3 

II 

o 

o 

bars 

The  WP  layer  has  been  given  some  strength  and  Yq  has  been  in- 
creased in  all  layers. 

3.  High  strength,  y = g(V,e):  This  calculation  duplicates 
calculation  (2)  with  the  cavity  equation  of  state  described  in 
Section  2.1. 

4.  High  sound  speed:  This  calculation  repeats  calcula- 
tion '(3)  with  higher  compressional  wave  velocities  in  each 
layer.  The  new  sound  speeds  correspond  closer  to  logs  from 
hole  U8al0. ^ To  reach  these  higher  velocities,  each  layer 
is  given  a larger  bulk  and  shear  modulus.  However,  its 
Poisson's  ratio  remains  the  same  as  for  the  corresponding 
layer  of  calculation  (3) . The  higher  sound  speeds  are 

Layer  1 CQ  = 2.03  km/sec  (6700  ft/sec) 

Layer  2 CQ  = 2.41  km/sec  (7900  ft/sec) 

Layer  3 CQ  = 2.0  km/sec  (6600  ft/sec) 

5.  Constant  Poisson's  ratio:  This  calculation  is  identi- 
cal to  calculation  (3)  except  that  the  shear  modulus  for  each 
calculational  zone  is  allowed  to  vary  during  porous  crushup  in 
order  to  maintain  Poisson's  ratio  constant  for  that  zone. 

The  results  of  this  study  are  summarized  in  Table  4. 
Comparing  calculations  (1)  and  (la) , clearly  the  thickness  of 
layer  2 is  not  a significant  factor  in  magnitude  of  the  peak 
residual  stresses  (note  that  all  residual  stresses  are  relative 
to  an  overburden  pressure  PQ  of  60  bars) , final  cavity  radius 
or  pressure,  or  RDP . The  RDP  (or  steady  state  value  of  the 
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reduced  displacement  potential)  is  defined  as  the  final  dis- 
placement multiplied  by  the  square  of  the  initial  radial 
distance  from  the  WP.  It  is  a constant  with  position  in  the 
far  (elastic)  field  and  is  a qualitative  measure  of  the  tele- 
seismic  coupling.  It  will  be  discussed  later  when  Baneberry 
results  are  compared  with  calculations  of  events  in  other  NTS 
media. 

Calculation  (2)  shows  that,  for  higher  material  strength, 
the  residual  stresses  are  increased  dramatically.  This  is  in 
agreement  with  the  results  of  Rimer, ^ which  indicate  that  in- 
creasing Yq  is  most  effective  in  increasing  residual  stress 
magnitudes.  Cavity  pressure  is  relatively  unchanged. 

When  the  cavity  equation  of  state  is  altered  as  in  cal- 
culation (3),  a larger  cavity  pressure  is  noted,  together  with 
a smaller  cavity  radius.  This  larger  cavity  pressure  is  still 
somewhat  smaller  than  overburden.  Modifications  in  compres- 
sional  wave  velocity  (calculation  4)  or  in  elastic  behavior 
(calculation  (5))  have  little  effect. 

This  sensitivity  study  indicates  that  the  residual 
stresses  are  primarily  affected  by  changes  in  material  strength 
while  cavity  pressure  is  modified  by  changes  in  the  cavity 
equation  of  state.  The  basic  results  of  Reference  4 remain 
unchanged.  The  residual  stresses,  while  larger  for  the  high 
strength  calculations,  are  still  much  smaller  than  the  results 
for  Yucca  Flat  tuff  and  for  Rainier  Mesa  tuff  presented  in 
Reference  4. 
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4.  THE  EFFECT  OF  CLAY  AND  WATER  ON  RESIDUAL  STRESSES 


Spherical,  one-dimensional  SKIPPER  calculations  for  a 
working  point  medium  having  extremely  high  clay  and  water 
content  (the  Baneberry  event)  are  discussed  in  this  section 
and  compared  with  calculations  for  Yucca  Flat  dry  tuff  and 
for  the  saturated  tunnel  tuffs  of  Rainier  Mesa.  All  calcula- 
tions use  the  same  constitutive  modeling  and  the  same  device 
yield  so  that  differences  in  the  results  are  a function  only 
of  the  material  properties  chosen  (constitutive  models  are 
discussed  in  detail  in  Section  2 and  material  properties  for 
the  calculations  are  given  in  Table  2) . Table  5 shows  the 
important  results  of  these  calculations. 

The  events  considered  for  Rainier  Mesa  are  Dido  Queen 
in  a weaker  tunnel  tuff  environment  and  Mighty  Epic  in  a 
stronger,  atypical  tunnel  tuff.  A second  Dido  Queen  calcula- 
tion was  made  for  comparison  using  the  Baneberry  overburden 
pressure  of  60  bars  which  barely  perturbed  the  results.  The 
Yucca  Flat  dry  tuff  calculations  represent  average  properties 
for  that  area.  These  two  Yucca  Flat  calculations  differ  only 
in  the  failure  envelope  used  as  do  the  two  Baneberry  calcula- 
tions presented.  All  Baneberry  calculations  assumed  2 layers; 
one  having  greater  than  50  percent  montmorillonite  clay,  and 
the  other  less  than  50  percent. 

Table  5 gives  the  maximum  values  of  the  residual 
stresses  in  the  radial  and  transverse  (hoop)  directions 
relative  to  the  scalar  overburden  pressure.  The  Baneberry 
calculations  showed  much  lower  residual  stresses  when  compared 
with  the  other  Yucca  Flat  calculations.  While  the  comparison 
with  the  weaker  Rainier  Mesa  calculations  (Dido  Queen)  was 
not  as  dramatic,  the  calculations  still  differ  considerably. 
The  most  important  material  properties  in  determining  the 
values  of  the  residual  stresses  are  the  magnitude  and  shape 
of  the  failure  surface.  Since  the  "strong"  Baneberry 
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Table 


(weak) 


calculation  uses  a failure  surface  determined  for  Diagonal 
Line  alluvium,  recognized  as  a very  strong,  atypical  alluvium, 
the  "weak"  calculation  is  probably  a more  believable  estimate 
of  the  residual  stresses  around  the  Baneberry  cavity. 

The  residual  stresses  must  prevent  the  cavity  pressure 
from  hydrofracturing  the  surrounding  rock.  Table  5 shows 
cavity  pressures  obtained  from  the  SKIPPER  calculations. 

These  cavity  pressures  are  the  least  reliable  part  of  the 
calculations  since  they  are  obtained  using  a simple  bubble 
model  described  in  Section  2 to  initiate  the  explosion. 

Cavity  radius  is  presented  for  each  calculation  at  the  same 
yield.  Again,  the  most  significant  influence  on  cavity  size 
is  the  failure  envelope,  higher  strength  giving  smaller 
cavities  and  vice  versa.  There  is  a general  trend  of  smaller 
cavities  indicating  higher  cavity  pressures. 

A primitive  measure  of  the  teleseismic  coupling  implied 
by  these  calculations  is  given  by  the  steady  state  value  of 
the  reduced  displacement  potential  (RDP)  shown  for  yields  cor- 
responding to  Baneberry.  Data  compiled  by  Savino^8^  of  S8 
show  that  both  for  regional  and  teleseismic  measurements, 
Baneberry  has  the  same  body  wave  magnitude  as  Yucca  Flat 
shots  of  4 to  5 times  the  Baneberry  yield.  Also,  for  the  same 
depth  of  burial,  Baneberry  has  4 or  5 times  the  surface  wave 

[19] 

amplitude.  The  data  based  on  Bache,  et  al.  therefore 

indicates  that  the  average  Yucca  Flat  calculation  has  the 
proper  RDP  relative  to  the  best  guess  "weak"  Baneberry  cal- 
culation. (RDP  varies  linearly  with  yield.)  Data  from 
Savino  also  shows  body  wave  magnitudes  from  Baneberry  com- 
parable to  typical  Rainier  Mesa  events.  This  tends  to  indi- 
cate that  either  Baneberry  calculation  would  be  satisfactory, 
with  the  "strong"  being  slightly  preferred. 

Other  Baneberry  calculations  were  made;  one  varying  the 
magnitude  of  the  yield  strength  in  layer  1 (the  high  clay 
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content  WP  layer) , and  another  the  location  of  the  interface 
between  layers  1 and  2.  These  calculations  showed  that  the 
strength  and  location  of  layer  1 had  no  influence  on  the  peak 
residual  stress  v/hich  occurred  in  layer  2 for  all  calculations. 
The  strength  and  location  of  layer  1 did  influence  cavity  size 
and  therefore  cavity  pressure. 

Figures  2-7  show  the  complete  residual  stress  fields 
for  the  calculations  of  Table  5 (Dido  Queen  1 is  not  included 
since  it  is  almost  identical  to  Dido  Queen  2) . Note  that  the 
overburden  pressure  has  not  been  subtracted  out  from  these 
stresses.  It  is  observed  that  higher  strength  not  only  in- 
creases the  peak  residual  stress,  it  also  makes  the  peak  values 
occur  closer  to  the  WP.  (This  is  partially  a consequence  of 
smaller  cavity  size.)  For  lower  strengths,  the  stress  field 
is  spread  out  over  a considerably  greater  volume,  however. 

All  the  calculations  considered  the  medium  to  be  infinite  in 
extent  and,  so,  neglected  free  surface  effect.  However,  the 
free  surface  will  have  the  effect  of  reducing  the  magnitude 
of  these  residual  stresses.  For  Baneberry,  where  the  free 
surface  is  only  278  meters  from  the  working  point,  well  within 
the  calculated  residual  stress  field,  the  residual  stress 
field  would  be  greatly  affected.  Although  these  results  can- 
not be  considered  definitive  evidence,  they  are  instructive 
and  thought  provoking.  They  strongly  suggest  that  devices 
should  be  buried  more  deeply  in  a weaker  medium. 
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Figure  4.  Residual  stress  field  versus  position  for  Yucca 
Flat  1 (weak) . 
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Figure  5.  Residual  stress  field  versus  position  for  Yucca 
Flat  2 (strong) . 
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Figure  6.  Residual  stress  field  versus  position  for  Dido 
Queen. 
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5.  FIGURES  OF  MERIT 


It  is  difficult  to  quantify  what  is  important  to  con- 
tainment of  the  cavity  gases  when  comparing  the  calculations 
of  Table  5 (Section  4).  Some  possible  "figures  of  merit"  are 
discussed  in  this  section.  In  Reference  4,  and  in  Section  4 
of  this  report,  the  values  of  the  peak  residual  stresses 
(particularly  the  hoop  stresses) , relative  to  the  scalar 
overburden  pressure,  were  compared.  It  is  reasonable  to  use 
the  stresses  with  the  overburden  subtracted  out  as  a measure 
of  containment  since  the  philosophy  is  that  it  is  the  addi- 
tional stresses  over  the  overburden  which  forms  the  contain- 
ment membrane  to  prevent  hydrofracture  from  the  cavity.  There 
it  was  stated  that  containment  was  better  when  the  peak 
residual  stresses  were  higher. 

This  simple  figure  of  merit  does  not  give  the  complete 
picture,  however.  Figures  2-7  show  that  the  residual  stress 
fields  look  quite  different  for  strong  media  (where  the  peak 
stresses  are  high)  and  weak  media.  For  weak  media,  the  com- 
pressive hoop  stresses,  though  smaller  in  magnitude,  tend  to 
occur  over  a larger  radial  distance,  i.e.,  the  residual 
stress  membrane  is  thicker.  However,  we  believe  that  a thin- 
ner, but  significantly  higher  stress  membrane  is  far  more  ef- 
fective as  a containment  membrane,  so  that  we  will  continue 
to  use  the  peak  residual  stress  in  our  figure  of  merit. 

Of  course,  these  residual  stresses  must  resist  a cavity 
pressure.  Based  on  the  simple  mechanics  of  materials  problem 
of  a thin  walled  pressure  vessel,  the  ratio  of  cavity  pressure 
to  peak  transverse  residual  stress  would  appear  to  be  a 
reasonable  figure  of  merit.  However,  the  cavity  pressure  com- 
puted using  SKIPPER  is  not  a good  measure  of  the  real  cavity 
pressure.  It  takes  no  account  of  water  getting  into  the 
cavity  from  the  surrounding  rock  or  of  additional  thermal 
equilibration  expected  in  the  rock  outside  the  calculated 
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cavity.  A procedure  developed  by  Peterson  and  Lie^20^ 
tain  cavity  pressures  and  temperatures  for  cavity  leak  cal- 
culations was  modified  here  and  used  to  calculate  a better 
value  of  cavity  pressure  from  the  output  of  the  SKIPPER  cal- 
culations. 

Starting  from  the  final  state  obtained  from  the 
SKIPPER  output  (cavity  radius,  cavity  energy,  internal  energy 
and  mass  of  adjacent  zones,  etc.)  and  beginning  at  the  cavity 
boundary,  small  amounts  of  mass  of  rock  (and  its  accompanying 
internal  energy)  are  equilibrated  with  the  melted  cavity 
material  which  is  at  a specific  internal  energy  greater  than 
the  melt  energy  of  the  rock.  This  calculational  procedure 
approximates  the  physical  processes  of  spall  and  equilibration. 
The  process  is  stopped  when  sufficient  mass  has  been  added 
so  that  the  equilibrated  "cavity"  is  at  the  melt  energy. 
Finally,  the  volume  of  steam  in  the  new  cavity  is  computed 
based  on  the  water  content  of  the  undisturbed  rock,  and  a 
final  cavity  pressure  is  calculated  (for  the  steam) . 

One  question  which  arises  involves  how  much  water  is 
in  the  undisturbed  rock.  The  water  content  by  weight,  f , 
given  in  Table  2 of  Section  2,  is  defined  as  the  amount  of 
water  that  can  be  removed  at  105°C.  This  does  not  include 
the  bound  water  in  the  rock.  For  Rainier  Mesa  tuff,  this 
amounts  to  approximately  7 percent  extra  water  by  weight. ^2(^ 
The  cavity  pressure  was  computed  for  each  of  the  calculations 
shown  in  Table  5,  using  fw  alone  and  also  including  the  ad- 
ditional bound  water  (assumed  to  be  7 percent  for  all  media) . 
Table  6 presents  the  results  of  these  computations.  Note  that 
the  effect  of  equilibrating  the  cavity  with  the  outside  mate- 
rial is  to  reduce  the  cavity  pressure  for  all  cases.  When  the 
bound  water  is  included  in  the  computation  of  cavity  pressure, 
less  of  a change  results,  as  would  be  expected.  (The  mass  of 
rock  material  reduces  the  pressure  while  the  water  tends  to 
increase  it.) 
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Table  6 

Figures  of  Merit  for  Equilibrated  and  Unequilibrated  Cavities 
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Here,  a is  the  peak  transverse  residual  stress  and 
the  figures  of  merit  chosen  are  ratios  of  the  calculated 
cavity  pressures,  P^,  P£  and  P^  to  ot>  Let  us  define  a higher 
ratio  of  cavity  pressure  to  a as  a greater  relative  contain- 
ment risk.  For  both  the  figures  of  merit  in  which  cavity 
pressure  from  an  equilibrated  cavity  was  used,  both  Baneberry 
calculations  are  clearly  less  conservative  and  therefore 
relatively  greater  containment  risks.  Using  the  unequilibrated 
cavity  pressure,  the  "weak",  best  guess,  Baneberry  2 calcula- 
tion is  the  worst  risk.  However,  the  "strong"  Baneberry  1 
calculation  lies  within  the  data  for  the  other  events. 

Summarizing,  for  all  figures  of  merit  investigated, 
the  "best  guess"  Baneberry  2 calculation  is  clearly  the  worst 
risk.  Only  when  the  unequilibrated  cavity  pressure  (a  hardly 
justifiable  choice)  is  used,  does  even  the  strong  upper  limit 
Baneberry  1 calculation  appear  as  safe  as  for  other  events. 
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THE  INFLUENCE  OF  CLAY  AND  WATER  IN  ROCKS 
ON  CONTAINMENT  OF  UNDERGROUND  NUCLEAR  EXPLOSIONS 

1.  INTRODUCTION  AND  SUMMARY 

This  report  discusses  the  effects  that  large  amounts  of 
clay  and  water  in  the  surrounding  rock  may  have  on  the  contain- 
ment of  hot  cavity  gases  due  to  an  underground  nuclear  event. 

The  containment  failure  of  the  Baneberry  event,  detonated  in 
Yucca  Flat  in  December  of  1970,  seems  to  be  an  ideal  vehicle 
for  this  study.  The  Baneberry  site  has  been  characterized  by 
the  USGS^  as  having  unusually  high  concentrations  of  mont- 
morillonite  clay  (65  percent)  around  the  working  point.  The 
water  content  of  the  working  point  material  (an  altered  tuff) 
is  related  to  the  clay  content  and  has  been  estimated  by  both 
the  USGS  ^ and  by  L.  Ramspott^^  to  be  as  much  as  25  percent 

by  weight.  Material  properties  data  for  Rainier  Mesa  tuff, 

T31 

compiled  by  S.  Butters  of  TerraTek,  Inc.  indicate  that  high 
water  content  in  laboratory  samples  correlates  with  low  shear 
strength.  While  there  is  little  hard  data  on  the  effect  of 
clay  on  material  strength,  it  seems  clear  that  the  presence 
of  clay  can  only  weaken  the  host  material. 

A scenario  for  the  containment  failure  of  Baneberry  has 
. r 4 1 

been  presented  in  a report  by  N.  Rimer  which  discusses  one- 
dimensional ground  motion  calculations  using  the  SKIPPER  code 
of  the  Baneberry  event  and  of  typical  Rainier  Mesa  and  Yucca 
Flat  events  which  were  contained.  Material  properties  for 
Baneberry  were  based  on  site  investigations  and  analyses  by 
USGS^  and  L.  Ramspott^  which  led  to  the  assumption  of  three 
spherical  layers  about  the  Baneberry  working  point.  The  working 
point  layer  of  altered  tuff  was  assumed  to  have  negligible  shear 
strength  due  to  the  very  high  clay  and  water  content.  The  second 
and  third  layers  for  the  calculations  were  both  alluvium  (the 
third  layer  was  of  infinite  extent) . Since  the  clay  content  of 
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the  second  layer  had  been  estimated  at  between  20  and  50  per- 
cent, it  was  assumed  to  have  half  the  shear  strength  of  the 
third  layer. 

The  major  difference  between  the  Baneberry  calculation 
and  the  Rainier  Mesa  and  Yucca  Flat  results  appeared  to  be  in 
the  magnitude  of  the  compressive  residual  stress  fields  due  to 
nonuniform  plastic  loading  and  subsequent  unloading  of  the  rock 
around  the  explosion-produced  cavities.  Residual  stress  fields 
are  seen  after  the  cavity  rebounds  in  all  ground  motion  calcu- 
lations at  S3.  Parameter  studies  by  N.  Rimer ^ have  indicated 
that  the  peak  transverse  residual  stresses  depend  very  strongly 
on  the  shear  strength  of  the  rock.  Therefore,  it  was  not  sur- 
prising that  the  Baneberry  calculation,  modeled  with  relatively 
low  strength,  gave  significantly  lower  residual  stresses. 

Based  on  these  results,  the  following  scenario  was  pre- 
sented as  possibly  leading  to  the  venting  of  cavity  gases  to 
the  surface  at  late  times.  The  low  strength  due  to  high  clay 
and  water  content  led  to  low  residual  stress  fields.  In  the 
absence  of  a significant  residual  stress  field,  the  hot  cavity 
gases  could  cause  tensile  hydrofracturing.  The  fracturing, 
if  sufficiently  great,  could  lead  to  venting.  The  low  material 
strength  was  considered  to  be  the  primary  reason  for  the  con- 
tainment failure  of  Baneberry. 

Although  the  original  report  of  this  work  was  written  in 

the  summer  of  1975,  it  received  only  limited  distribution  to 

interested  ERDA  representatives  at  the  time.  Because  of  the 

general  sensitivity  of  the  subject  matter,  we  were  asked  to 

% 

withhold  further  distribution  until  additional  calculations  were 
completed  which  would  determine  the  sensitivity  of  the  results 
to  the  modeling  assumptions.  A set  of  suggested  multi-layer 
calculations  was  completed  and  presented  at  an  Earth  Motion 
Calculators  meeting  held  in  La  Jolla  on  May  4,  1976.  Since  the 
multi-dimensional  calculations  requested  in  Ref.  4 have  recently 


been  done  and  presented  by  LLL  at  the  66th  CEP  meeting,  the 
report  was  released  in  January  of  1977. 

In  the  intervening  period,  considerable  work  has  been 
done  to  clarify  questions  raised  by  that  report.  The  present 
work  includes  a discussion  of  constitutive  modeling  and 
material  properties  data  improvements  (see  Section  2)  used  in 
later  calculations  of  both  Baneberry  and  the  comparison  events 
in  Rainier  Mesa  and  Yucca  Flat.  It  was  felt  that  the  com- 
parisons in  Ref.  4 might  not  be  completely  valid  since  there 
were  some  modeling  differences  in  the  calculations  being  com- 
pared (the  Rainier  Mesa  and  Yucca  Flat  calculations,  as  des- 
cribed in  Ref.  4,  were  not  made  specifically  as  part  of  this 
program.)  Careful  consideration  has  shown  that  those  comparisons 
are  qualitatively  valid. 

Section  3 describes  the  results  of  calculations  presen- 
ted at  the  Earth  Motion  Calculators  meeting  of  May  4,  1976. 

Based  on  suggestions  from  ERDA  representatives,  the  sensitivity 
of  the  results  to  some  of  the  modeling  assumptions  used  in  the 
Baneberry  calculation  was  investigated.  It  was  found  that  the 
basic  results  of  Ref.  4 were  not  sensitive  to  any  parameter 
other  than  material  strength. 

In  Section  4,  calculations  are  presented  for  a weak  and 
a strong  Rainier  Mesa  saturated  tuff,  a weak  and  a strong 
Yucca  Flat  dry  tuff,  as  well  as  a weak,  best-guess  Baneberry 
calculation  and  an  upper-limit,  relatively-high-strength  Bane- 
berry calculation.  These  calculations  show  Baneberry  to  have 
a much  lower  residual  stress  field  than  the  other  NTS  events 
due  to  the  lower  strength  of  the  surrounding  alluvium. 

Section  5 discusses  "figures  of  merit,"  measures  by 
which  one  can  tell  whether  a particular  event  is  a greater 
containment  risk  than  other  previously-detonated  successful 
events.  The  figure  of  merit  chosen  was  the  ratio  of  the  cavity 


pressure  to  the  peak  transverse  residual  stress.  For  this 
comparison,  a sophisticated  determination  of  the  cavity 
pressure  was  used  as  described  below.  Using  the  resulting 
figure  of  merit,  Baneberry  was  clearly  the  greatest  contain- 
ment risk  of  the  events  examined. 
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2.  CONSTITUTIVE  MODELING  ^ND  MATERIAL  PROPERTIES 


In  Section  2.1,  the  constitutive  models  are  discussed 
briefly  with  emphasis  on  the  modeling  improvements  between 
Ref.  4 and  the  calculations  presented  here.  Section  2.2  gives 
the  material  properties  used  for  the  calculations  of  Ref.  4 
and  of  Section  3 of  this  report.  The  modeling  and  material 
properties  used  for  the  results  presented  in  Section  4 are 
discussed  in  Section  2.3. 


2.1  CONSTITUTIVE  MODELING 

Descriptions  are  given  here  of  the  constitutive  models 
used  in  the  SKIPPER  calculations.  A more  complete  discussion 
may  be  found  in  Cherry,  et.  al . ^ 


Equations  of  State 

The  Baneberry  calculations  described  in  Ref.  4 used  an 
equation  of  state  of  the  form  developed  by  Tillotson^  and 
fit  at  S3  to  available  shock  data  for  a number  of  geological 
materials.  For  compressed  states  (p  > p ) and  for  cold  expan- 
ded states  (p  < p and  e < e ) , the  pressure  is  given  by 
o s 


P 

s 


a + 


ep  + Ap  + Bu2 


For  expanded  states  (p  < p ) where  e > e^  , the  pressure  is 
given  by 
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lit 


The  phase  transition  from  liquid  to  vapor  (p  < pQ  and 

e < e < e')  is  approximated  as 
s s 


p 


es) 


p + 
v 


(es 


Ps]  • 


Here 


p = mass  density 
n - p/pq 
w = n + l 

e = specific  internal  energy  as  the  material  is 
brought  to  vaporization  temperature 

e'  = additional  specific  internal  energy  required 
to  change  the  material  from  a liquid  to  a 
vapor  state 


All  the . above  formulas  assume  that  no  air-filled  voids  are 
present  in  the  rock.  The  constants,  a,  b,  and  eQ  govern  the 
energy  dependence  of  the  material,  allowing  it  to  behave  as 
an  ideal,  y-law  9as  at  high  energies  (a  ^ = y-1)  » and  reducing 
the  effective  y at  lower  energies.  The  density  dependence  is 
governed  by  a quadratic  polynomial  with  constants  A and  B 
which  may  be  fit  to  shock  data.  The  Tillotson  equation  of 
state  does  not  accurately  model  the  expansion  from  rock  states 
above  approximately  10  kbars  where  the  expansion  of  steam  is 
important. 


The  Tabular  Array  of  Mixtures  Equation  of  State  (TAMEOS) 
described  in  Cherry^ ^ has  been  used  for  all  tuff  calculations 
discussed  here.  TAMEOS  mixes  the  proper  percentage  of  water 
with  grain  density  rock,  assuming  pressure  equilibrium  between 
the  rock  and  the  water,  to  give  the  pressure  response  for  the 
fully  crushed  tuff  mixture.  This  gives  a far  better  description 
of  the  subsequent  expansion  of  the  shocked  water-rock  mixture. 
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Porous  Crushup 

The  equations  of  state  discussed  above  describe  the 
pressure  response  for  fully  saturated  rocks.  Air-filled 
porosity  <pQ  is  included  through  ,the  S3  porous  crushup  model 
(P-a  model)  ^ in  which  the  pressure  of  the  porous  material 
is  described  by 


where 

V = specific  volume  of  porous  material 

e = specific  internal  energy  of  porous  material 
P = pressure  obtained  from  equation  of  state 
a = distension  ratio  defined  by  V/V  > 1 

V = specific  value  of  the  material  with  zero  air- 

filled  voids 

The  distension  ratio  is  required  to  decrease  from  an 
initial  value 


at  zero  pressure  down  to  1.0  as  the  pressure  increases  to  Pc» 
the  crush  pressure  or  pressure  limit  at  which  all  air-filled 
porosity  is  assumed  irreversibly  removed.  Any  porosity  lost 
during  compression  below  the  elastic  pressure  Pg  is  recover- 
able upon  expansion.  Any  porous  crushup  occuring  between  Pg 

and  P„  is  irreversible, 
c 

Elastic  Properties  and  Material  Strength 

The  ambient  compressional  wave  velocity  CQ  is  related 
to  the  bulk  modulus  Kq  and  the  shear  modulus  G by  the 
relationship 
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2 4 

p C = K + 2.  G 
o o 3 


Here,  the  shear  modulus  is  assumed  constant.  Since  the  bulk 
modulus  K varies  with  pressure,  Poisson's  ratio  o is  not 
being  held  constant.  One  calculation  was  made  with  constant 
a throughout  the  crushup.  This  requires  a constant  ratio  of 
K/G,  i.e.,  a varying  G. 

The  deviatoric  components  of  the  stress  tensor  are 


related  to  the  deviatoric  strain  rates 


S.  . = 2 G e.  . 
13  13 


by  Hooke's  law 


The  material  strength  model  requires  that  the  principal  stress 
be  within  the  von  Mises  yield  surface.  For  spherical  geometry, 
this  reduces  to  the  condition  that  the  magnitude  of  the  radial 
deviatoric  stress  not  exceed  2/3  the  yield  strength  (maximum 
stress  difference) . Behavior  of  a failed  element  is  governed 
by  the  non-associated  flow  rule. 

For  the  calculations  of  Ref.  4,  the  yield  strength  was 
given  by 


* = k + ^ (2  - yi  i1  - §-) p < v 

Y = (y  + Y ) (l  - — ) P > P , 

\ ° »/  \ em/  “ m 


e < e 


P > P , e < e 
— m m 


Y = 0 


e > e 


where 


Y , Y , 
o m 


are  constants  for  a given  material. 


The  more  recent  calculations  use  an  improved  form  for 


the  yield  strength  which  has  been  shown  by  Cherry  and  Peterson 
to  be  empirically  superior  to  a simple  pressure  dependence. 
Here,  the  pressure  P in  the  above  equations  is  replaced  by 
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P which  is  a function  of  P and  the  deviatoric  stresses 
given  by 


P = P 


1 |S11S22S33 | 


1/3 


A description  of  how  shear  failure  is  treated  in  this  model  is 
given  in  Cherry,  et^.  al . ^ 

Tension  failure  is  allowed  to  occur  in  an  element  if  a 

principal  stress  becomes  tensile.  The  tension  failure  model 

r 91 

proposed  by  Maenchen  and  Sack  and  described  in  detail  by 
Cherry ^ is  then  applied.  An  inelastic  strain  which  is  just 
sufficient  to  zero  the  tensile  stress  is  introduced  normal  to 
the  crack.  This  strain  increases  or  decreases  as  the  crack 
opens  or  closes. 


Cavity  Equation  of  State 

The  cavity  source  region  was  chosen  to  be  initially 
large  enough  to  vaporize  70  metric  tons  of  rock  for  each  kilo- 
ton  of  device  yield.  (The  rock  is  assumed  to  be  at  undisturbed 
density.)  The  cavity  pressure  for  all  calculations  is  given  by 

P = (y-1)  e/V 

For  the  Rainier  Mesa  and  Yucca  Flat  calculations  des- 
cribed in  Ref.  4,  y was  set  equal  to  1.5.  The  Baneberry  cal- 
culation of  Ref.  4 assumed  that  y was  a function  of  specific 
volume  given  by 

y = 1.03  + 0.9//V 

This  expression  for  y was  a preliminary  result  obtained  for 
tuff  modeled  using  the  chemical  equilibrium  CHEST  code  of 
D.  Laird and  has  been  modified  for  later  calculations.  It 
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will  be  shown  later  that  the  peak  residual  stresses  are  not 
sensitive  to  the  choice  of  cavity  y.  However,  depending  upon 
the  cavity  model  used,  peak  velocities  at  a given  location  may 
differ  by  20  percent. 

More  recent  calculations  used  the  following  expression 

for  y. 


y ' = 1.085 

y'  = 1.2085  - 0.065  V 
y'  = 1.9375  - 0.875  V 
y = y'  - 0.032  + 3.2  x 10-13  e 

y = y' 


V > 1.9 

1.9  > V >_  0.9 

V < 0.9 

e > 1011  ergs/gm 
e < 1011 


It  should  be  emphasized  that  the  above  expressions  for 
y are  only  very  approximate  fits  to  the  complicated  tabular 
equation  of  state  of  Laird.  The  fit  given  above  was  derived 
from  data  from  a table  developed  for  one  particular  water 
content  (24  percent  by  weight,  including  approximately  7 per- 
cent bound  water) . Since  the  creation  of  a CHEST  table  is  a 
complicated  procedure,  the  same  expression  for  y was  used  in 
calculations  for  Baneberry,  Rainier  Mesa  and  Yucca  Flat  tuffs, 
i.e.,  over  a wide  range  of  water  contents. 


2.2  MATERIAL  PROPERTIES  FOR  THE  PARAMETER  STUDY 

Table  1 gives  the  material  properties  data  used  in  the 
original  calculations  of  Ref.  4.  For  Baneberry,  Material  1 
describes  a saturated,  high  clay  content  altered  tuff  extend- 
ing radially  15.24  meters  (50  feet)  from  the  working  point. 
Since  the  clay  content  of  this  layer  was  greater  than  50  per- 
cent, the  material  strength  was  considered  negligible.  The 
Tillotson  equation  of  state,  using  constants  approximating  the 
data  of  Stephens,  et.  al^.  was  used  to  model  the  pressure 

response  of  the  material. 
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Materials  2 and  3 are  assumed  to  be  alluvium  modeled 


using  a Tillotson  equation  of  state.  Material  2,  with  greater 
than  20  percent  montmorillonite , is  assumed  to  have  half  the 
shear  strength  of  Material  3,  which  is  a weak  alluvium. 

Material  2 was  assumed  to  be  36.58  meters  thick  (120  feet) 
and  to  be  surrounded  by  Material  3.  Free  surface  effects  and 
a radially  varying  gravitational  field  were  neglected  in  these 
calculations. 

[1  21 

The  measured  data  ' upon  which  the  modeling  was  based 
was  often  conflicting  since  the  Baneberry  geology  is  unusually 
complicated  and  since  the  data  were  accumulated  from  holes 
drilled  both  pre  and  postshot  and  in  different  locations  rela- 
tive to  the  WP.  The  log  data  indicated  that  Material  2 had  a 
compressional  wave  velocity  20  to  30  percent  greater  than 
Materials  1 or  3.  This  was  included  in  the  modeling.  Also, 
due  to  the  higher  clay  content,  Material  2 had  only  half  the 
air-filled  void  content  of  Material  3.  Also  included  in  Table  1 
is  material  properties  data  used  for  calculations  for  Yucca 
Flat  dry  tuff  and  for  Rainier  Mesa  saturated  tuff.  These  cal- 
culations were  not  made  specifically  for  the  Baneberry  study 
but  were  used  for  comparison  with  the  Baneberry  results.  As  a 
result,  differences  exist  in  the  constitutive  modeling  for  the 
three  calculations.  For  example,  both  Yucca  Flat  and  Rainier 
Mesa  calculations  use  a TAMEOS  equation  of  state  rather  than  a 
Tillotson  equation  of  state. 

A parameter  study  was  made  to  examine  the  sensitivity  of 
the  results  of  Ref.  4 to  the  modeling  assumptions  used  in  the 
Baneberry  calculation.  The  starting  point  for  this  study  is 
the  material  properties  data  of  Table  1.  The  calculations  of 
Ref.  4 all  used  the  pressure  dependent  failure  surface  described 
in  Section  2.1.  However,  the  parameter  study  used  the  empiri- 
cally superior  P dependent  failure  surface.  The  parameter  study 
also  used  the  newer  functional  dependence  of  y on  specific  volume 


■ 

l 

, 


and  energy  which  was  described  in  Section  2.1  for  the  cavity 
equation  of  state . 

2.3  MATERIAL  PROPERTIES  FOR  CALCULATIONS  OF  SECTION  4 

Using  the  constitutive  modeling  described  in  Section  2.1, 
a series  of  calculations  was  made  for  Rainier  Mesa  (Dido  Queen 
and  Mighty  Epic  events) , for  Yucca  Flat  dry  tuff  (a  calcula- 
tion using  average  properties)  and  for  Baneberry.  Table  2 
gives  the  properties  used  for  these  calculations.  All  calcu- 
lations use  the  TAMEOS  equation  of  state  for  the  working  point 
material , a cavity  equation  of  state  where  y is  a function 
of  specific  volume  and  energy,  and  the  P dependence  for  the 
failure  surface.  This  section  discusses  the  properties  chosen 
for  each  event. 

The  material  properties  for  the  Dido  Queen  event,  den- 
sity, grain  density,  water  content  and  compressional  wave 
velocity  were  obtained  from  the  CEP  document  for  that  event. 

The  failure  envelope  was  modeled  using  Figure  24  of  Butters, 

[ 12] 

et.  al_.  while  the  crush  curve  was  obtained  from  Figure  19 

of  the  same  report.  Dido  Queen  was  chosen  for  this  study 
because  it  was  located  in  tunnel  e of  the  Rainier  Mesa  complex, 
a tunnel  that  triaxial  loading  tests  indicate  contains  rela- 
tively weaker  tuffs. 

The  Mighty  Epic  event,  located  in  tunnel  n,  was  chosen 
as  an  example  of  a stronger  Rainier  Mesa  tuff.  The  choice  of 
material  properties  for  the  Mighty  Epic  event  is  discussed  in 
detail  in  Rimer,  et^.  al . 

The  Yucca  Flat  material  properties  were  based  on  average 

, r 14 1 

properties  for  a Yucca  Flat  dry  tuff  summarized  by  Ramspott. 

Unfortunately,  there  is  no  data  on  the  strength  of  this  mate- 
rial. TerraTek  laboratory  data  for  Rainier  Mesa  tunnel  bed 
tuffs  compiled  by  S.  Butters ^ shows  a relationship  between 
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shear  strength  and  water  content.  It  should  be  emphasized 
that  the  data  of  Ref.  3 is  stress  difference  in  uniaxial  strain 
at  4 kbars  confining  pressure,  not  the  yield  strength.  However, 
it  should  give  a relative  indicator  of  the  strength  of  the 
materials. 

F.  App^15^  has  drawn  several  straight  lines  through  a 
plot  of  strength  (stress  difference  of  4 kbars)  v£  water  con- 
tent (by  weight)  for  these  data  which  indicate  that  the  shear 
strength  is  inversely  proportional  to  the  water  content. 

These  lines  shown  in  Figure  1 imply  that  for  Baneberry  water 
contents,  the  strength  is  quite  small.  For  the  lower  water 
content  of  Yucca  Flat  tuff  above  the  water  table,  the  plot 
gave  a strength  value  of  0.63  kb.  The  yield  strength  of 
0.63  for  Yucca  Flat  tuff  water  content  was  apportioned  between 
Yq  and  Y^  in  the  ratio  suggested  by  the  Dido  Queen  data.  This 

leaves  the  magnitude  of  P yet  to  be  determined.  For  this 

m 

reason,  two  calculations  were  made  for  Yucca  Flat,  one  for  a 

P„  of  3.0  kbar,  as  in  Dido  Queen,  and  another  for  a P of  1.0 
m m 

kbar,  as  in  Mighty  Epic.  The  parameter  Pm  is  as  important  to 
the  residual  stress  field  as  the  magnitude  of  the  failure 
strength  itself.  A low  value  of  Pm  indicates  a steeper  slope 
of  the  failure  surface  and  therefore  a higher  strength  at  low 
pressures. 

A Yucca  Flat  calculation  made  with  the  data  of  Table  2 
using  the  TAMEOS  equation  of  state  was  compared  to  one  made 
using  the  Tillotson  equation  of  state.  Tillotson  constants 
A and  B for  this  calculation  were  chosen  to  match  the  zero 
pressure  bulk  modulus  used  in  the  TAMEOS  calculation  and  to 
match  its  zero  energy  curve  at  the  crush  pressure  P . The 
comparisons  are  given  in  Table  3. 
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Uniaxial  stress  difference  at  4 kbars  (range  and 
mean)  versus  water  content  (by  weight)  for  Terra  Tek 
laboratory  data. 


Table  3 


Comparison  of  results  using  TAMEOS  and  Tillotson 
equations  of  state  for  Yucca  Flat  dry  tuff 


Cavity 

Cavity 

Radial 

Transverse 

radius 

pressure 

stress 

stress 

(meters) 

(bars) 

(bars) 

(bars) 

TAMEOS 

31.70 

136.7 

174 

304 

Tillotson 

29.18 

183.9 

140 

249 

The 

calculations  showed  about 

a 10  percent  difference 

in  cavity 

radius  and 

close  to  20  percent  differences  in 

radial  and 

transverse 

peak  residual 

stresses. 

Cavity  pres- 

sures  were  even  further  apart.  On  the  basis  of  these  results, 
it  was  decided  to  apply  TAMEOS  to  the  Baneberry  working  point 
material.  Input  to  TAMEOS  is  the  water  content  (by  weight) , 
the  grain  density  of  the  tuff,  and  constants  defining  the 
equation  of  state  of  the  grain  density  tuff.  It  was  assumed 
that  the  equation  of  state  constants  applicable  to  Rainier  Mesa 
tuff  (having  much  lower  grain  densities)  may  be  used  for  the 
Baneberry  altered  tuff.  The  important  features  are  the  high 
water  content  and  the  high  grain  density  due  to  clay  content. 

Since  the  peak  residual  stresses  had  been  shown  to  be 

primarily  dependent  on  the  failure  strength  and  not  to  be  very 

dependent  upon  the  exact  position  of  the  layer  interface  or  on 

the  compressional  wave  velocity  (the  results  are  presented  in 

Section  3) , it  was  decided  to  model  the  Baneberry  site  using 

only  two  layers,  a working  point  layer  of  over  50  percent 

montmorillonite  extending  radially  35.66  meters  (117  feet)  from 

the  WP,  with  an  alluvium  layer  of  lower,  but  still  significant 

clay  content  around  it.  The  location  of  the  layer  was  based 

r 2 ] 

on  data  from  drill  hole  Ue8i  discussed  by  Ramspott.  The 

Tillotson  equation  of  state  was  used  for  the  second  layer  which 
is  not  expected  to  see  stresses  greater  than  5-10  kbars.  Both 
layers  were  modeled  using  data  from  hole  Ue8i. 
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For  layer  2,  data  on  the  failure  envelope  has  been 

found  for  Diagonal  Line  alluvium, considered  to  be  a rela- 

f 17 1 

tively  strong  alluvium  and  for  Merlin  alluvium,  considered 

to  be  a more  representative  alluvium.  The  strength  given  for 
layer  2 for  the  calculation  labeled  Baneberry  1 in  Table  2 
was  obtained  using  Diagonal  Line  alluvium  data  and  so  repre- 
sents an  upper  limit  to  the  strength  of  layer  2.  The  calcu- 
lation labeled  Baneberry  2 is  identical  with  Baneberry  1 
except  that  the  strength  of  layer  2 is  lower  as  shown  since 
it  was  modeled  using  Merlin  data. 

The  failure  strength  of  the  Baneberry  working  point 
layer  was  modeled  using  data  from  Stephens,  et.  al. 

Figure  7 of  Ref.  17  plots  shear  strength  V£  water  contents  for 
alluvium.  That  plot  indicates  that  for  water  contents  of  20 
percent  or  more,  alluvium  has  zero  shear  strength.  This  implies 
that  the  working  point  material  may  have  negligible  strength. 
Figure  1 indicates  negligible  strength  for  water  contents  of 
30  percent. 
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3.  SENSITIVITY  OF  BANEBERRY  RESULTS  TO  MODELING  ASSUMPTIONS 


The  results  of  a parameter  study  presented  at  the 
Earth  Motion  Calculators  meeting  held  in  La  Jolla  on  May  4, 
1976  are  discussed  in  this  section.  The  purpose  of  this 
study  was  to  answer  questions  raised  by  ERDA  representatives 
as  to  the  sensitivity  of  the  results  of  Reference  4 to  the 
modeling  assumptions  used  for  the  Baneberry  calculation.  The 
material  properties  used  for  the  3 material  Baneberry  base 
calculation  are  those  used  for  Reference  4 and  are  given  in 
Table  1 of  Section  2.2.  A second  look  at  the  modeling 
assumptions  indicated  that  the  thickness  of  layer  2 (mate- 
rial 2)  should  be  51.8  meters  (170  feet)  rather  than  36.58 
meters  (120  feet)  as  in  Reference  4.  Modeling  improvements 
were  incorporated  into  this  study  as  described  in  Section  2. 

n_  The  calculations  of  this  study  are  the  following: 

1.  The  base  calculation:  This  calculation  uses  the  data 

of  Table  1 but  assumes  that  layer  2 is  170  feet  thick.  An 
important  feature  is  that  the  cavity  gas  equation  of  state 
has  a constant  y of  1.5. 

The  yield  strengths  are 

Layer  1 (clay)  Y = 0 

Layer  2 (half  strength)  Yq  = 7.0  bars,  Y^  = 50  bars, 

Pm  = 700  bars 
m 

Layer  3 (full  strength)  Yq  = 7.0  bars,  Y^  = 100  bars, 

Pm  = 700  bars 
m 

la.  This  calculation  is  identical  with  the  base  calcula- 
tion except  that  layer  2 is  120  feet  thick.  It  essentially 
duplicates  the  Baneberry  calculation  of  Reference  4 with 
better  modeling  for  the  failure  surface  and  with  a constant 
y for  the  cavity  gas. 
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2.  High  strength,  constant  y:  This  calculation  is 

identical  with  calculation  (1)  except  for  material  strengths. 
The  new  strengths  are 


Layer 

1 

Y = 
o 

30 

bars, 

Ym  = 10 

bars. 

P = 700 
m 

bars 

Layer 

2 

Y = 
o 

30 

bars. 

Ym  “ 50 

bars. 

P = 700 
m 

bars 

Layer 

3 

Y = 
o 

30 

bars. 

Ym  = 100 

bars 

, Pm  = 700 

bars 

The  WP  layer  has  been  given  some  strength  and  Yq  has  been  in- 
creased in  all  layers. 

3.  High  strength,  y = g(V,e):  This  calculation  duplicates 
calculation  (2)  with  the  cavity  equation  of  state  described  in 
Section  2.1. 

4.  High  sound  speed:  This  calculation  repeats  calcula- 
tion (3)  with  higher  compressional  wave  velocities  in  each 
layer.  The  new  sound  speeds  correspond  closer  to  logs  from 
hole  U8al0. ^ To  reach  these  higher  velocities,  each  layer 
is  given  a larger  bulk  and  shear  modulus.  However,  its 
Poisson's  ratio  remains  the  same  as  for  the  corresponding 
layer  of  calculation  (3) . The  higher  sound  speeds  are 


Layer  1 

CQ  = 2.03  km/sec 

(6700  ft/sec) 

Layer  2 

C =2.41  km/sec 

G 

(7900  ft/sec) 

Layer  3 

CQ  = 2.0  km/sec 

(6600  ft/s ec) 

5. 

Constant 

Poisson's  ratio: 

This  calculation  is  identi 

cal  to  calculation  (3)  except  that  the  shear  modulus  for  each 
calculational  zone  is  allowed  to  vary  during  porous  crushup  in 
order  to  maintain  Poisson's  ratio  constant  for  that  zone. 

The  results  of  this  study  are  summarized  in  Table  4. 
Comparing  calculations  (1)  and  (la) , clearly  the  thickness  of 
layer  2 is  not  a significant  factor  in  magnitude  of  the  peak 
residual  stresses  (note  that  all  residual  stresses  are  relative 
to  an  overburden  pressure  PQ  of  60  bars) , final  cavity  radius 
or  pressure,  or  RDP . The  RDP  (or  steady  state  value  of  the 
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reduced  displacement  potential)  is  defined  as  the  final  dis- 
placement multiplied  by  the  square  of  the  initial  radial 
distance  from  the  WP.  It  is  a constant  with  position  in  the 
far  (elastic)  field  and  is  a qualitative  measure  of  the  tele- 
seismic  coupling.  It  will  be  discussed  later  when  Baneberry 
results  are  compared  with  calculations  of  events  in  other  NTS 
media. 

Calculation  (2)  shows  that,  for  higher  material  strength, 
the  residual  stresses  are  increased  dramatically.  This  is  in 
agreement  with  the  results  of  Rimer, ^ which  indicate  that  in- 
creasing Yq  is  most  effective  in  increasing  residual  stress 
magnitudes.  Cavity  pressure  is  relatively  unchanged. 

When  the  cavity  equation  of  state  is  altered  as  in  cal- 
culation (3),  a larger  cavity  pressure  is  noted,  together  with 
a smaller  cavity  radius.  This  larger  cavity  pressure  is  still 
somewhat  smaller  than  overburden.  Modifications  in  compres- 
sional  wave  velocity  (calculation  4)  or  in  elastic  behavior 
(calculation  (5))  have  little  effect. 

This  sensitivity  study  indicates  that  the  residual 
stresses  are  primarily  affected  by  changes  in  material  strength 
while  cavity  pressure  is  modified  by  changes  in  the  cavity 
equation  of  state.  The  basic  results  of  Reference  4 remain 
unchanged.  The  residual  stresses,  while  larger  for  the  high 
strength  calculations,  are  still  much  smaller  than  the  results 
for  Yucca  Flat  tuff  and  for  Rainier  Mesa  tuff  presented  in 
Reference  4. 


4.  THE  EFFECT  OF  CLAY  AND  WATER  ON  RESIDUAL  STRESSES 


Spherical,  one-dimensional  SKIPPER  calculations  for  a 
working  point  medium  having  extremely  high  clay  and  water 
content  (the  Baneberry  event)  are  discussed  in  this  section 
and  compared  with  calculations  for  Yucca  Flat  dry  tuff  and 
for  the  saturated  tunnel  tuffs  of  Rainier  Mesa.  All  calcula- 
tions use  the  same  constitutive  modeling  and  the  same  device 
yield  so  that  differences  in  the  results  are  a function  only 
of  the  material  properties  chosen  (constitutive  models  are 
discussed  in  detail  in  Section  2 and  material  properties  for 
the  calculations  are  given  in  Table  2) . Table  5 shows  the 
important  results  of  these  calculations. 

The  events  considered  for  Rainier  Mesa  are  Dido  Queen 
in  a weaker  tunnel  tuff  environment  and  Mighty  Epic  in  a 
stronger,  atypical  tunnel  tuff.  A second  Dido  Queen  calcula- 
tion was  made  for  comparison  using  the  Baneberry  overburden 
pressure  of  60  bars  which  barely  perturbed  the  results.  The 
Yucca  Flat  dry  tuff  calculations  represent  average  properties 
for  that  area.  These  two  Yucca  Flat  calculations  differ  only 
in  the  failure  envelope  used  as  do  the  two  Baneberry  calcula- 
tions presented.  All  Baneberry  calculations  assumed  2 layers; 
one  having  greater  than  50  percent  montmorillonite  clay,  and 
the  other  less  than  50  percent. 

Table  5 gives  the  maximum  values  of  the  residual 
stresses  in  the  radial  and  transverse  (hoop)  directions 
relative  to  the  scalar  overburden  pressure.  The  Baneberry 
calculations  showed  much  lower  residual  stresses  when  compared 
with  the  other  Yucca  Flat  calculations.  While  the  comparison 
with  the  weaker  Rainier  Mesa  calculations  (Dido  Queen)  was 
not  as  dramatic,  the  calculations  still  differ  considerably. 
The  most  important  material  properties  in  determining  the 
values  of  the  residual  stresses  are  the  magnitude  and  shape 
of  the  failure  surface.  Since  the  "strong"  Baneberry 
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Table 


(weak 


calculation  uses  a failure  surface  determined  for  Diagonal 
Line  alluvium,  recognized  as  a very  strong,  atypical  alluvium, 
the  "weak"  calculation  is  probably  a more  believable  estimate 
of  the  residual  stresses  around  the  Baneberry  cavity. 

The  residual  stresses  must  prevent  the  cavity  pressure 
from  hydrofracturing  the  surrounding  rock.  Table  5 shows 
cavity  pressures  obtained  from  the  SKIPPER  calculations. 

These  cavity  pressures  are  the  least  reliable  part  of  the 
calculations  since  they  are  obtained  using  a simple  bubble 
model  described  in  Section  2 to  initiate  the  explosion. 

Cavity  radius  is  presented  for  each  calculation  at  the  same 
yield.  Again,  the  most  significant  influence  on  cavity  size 
is  the  failure  envelope,  higher  strength  giving  smaller 
cavities  and  vice  versa.  There  is  a general  trend  of  smaller 
cavities  indicating  higher  cavity  pressures. 

A primitive  measure  of  the  teleseismic  coupling  implied 
by  these  calculations  is  given  by  the  steady  state  value  of 
the  reduced  displacement  potential  (RDP)  shown  for  yields  cor- 
responding to  Baneberry.  Data  compiled  by  Savino^8^  of  S8 
show  that  both  for  regional  and  teleseismic  measurements, 
Baneberry  has  the  same  body  wave  magnitude  as  Yucca  Flat 
shots  of  4 to  5 times  the  Baneberry  yield.  Also,  for  the  same 

depth  of  burial,  Baneberry  has  4 or  5 times  the  surface  wave 

[191 

amplitude.  The  data  based  on  Bache,  et  al.  therefore 

indicates  that  the  average  Yucca  Flat  calculation  has  the 
proper  RDP  relative  to  the  best  guess  "weak"  Baneberry  cal- 
culation. (RDP  varies  linearly  with  yield.)  Data  from 
Savino  also  shows  body  wave  magnitudes  from  Baneberry  com- 
parable to  typical  Rainier  Mesa  events.  This  tends  to  indi- 
cate that  either  Baneberry  calculation  would  be  satisfactory, 
with  the  "strong"  being  slightly  preferred. 

Other  Baneberry  calculations  were  made;  one  varying  the 
magnitude  of  the  yield  strength  in  layer  1 (the  high  clay 
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content  WP  layer) , and  another  the  location  of  the  interface 
between  layers  1 and  2.  These  calculations  showed  that  the 
strength  and  location  of  layer  1 had  no  influence  on  the  peak 
residual  stress  which  occurred  in  layer  2 for  all  calculations. 
The  strength  and  location  of  layer  1 did  influence  cavity  size 
and  therefore  cavity  pressure. 

Figures  2-7  show  the  complete  residual  stress  fields 
for  the  calculations  of  Table  5 (Dido  Queen  1 is  not  included 
since  it  is  almost  identical  to  Dido  Queen  2) . Note  that  the 
overburden  pressure  has  not  been  subtracted  out  from  these 
stresses.  It  is  observed  that  higher  strength  not  only  in- 
creases the  peak  residual  stress,  it  also  makes  the  peak  values 
occur  closer  to  the  WP.  (This  is  partially  a consequence  of 
smaller  cavity  size.)  For  lower  strengths,  the  stress  field 
is  spread  out  over  a considerably  greater  volume,  however. 

All  the  calculations  considered  the  medium  to  be  infinite  in 
extent  and,  so,  neglected  free  surface  effect.  However,  the 
free  surface  will  have  the  effect  of  reducing  the  magnitude 
of  these  residual  stresses.  For  Baneberry,  where  the  free 
surface  is  only  278  meters  from  the  working  point,  well  within 
the  calculated  residual  stress  field,  the  residual  stress 
field  would  be  greatly  affected.  Although  these  results  can- 
not be  considered  definitive  evidence,  they  are  instructive 
and  thought  provoking.  They  strongly  suggest  that  devices 
should  be  buried  more  deeply  in  a weaker  medium. 
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4.  Residual  stress  field  versus  position  for  Yucca 
Flat  1 (weak) . 
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Figure  5.  Residual  stress  field  versus  position  for  Yucca 
Flat  2 (strong) . 
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5.  FIGURES  OF  MERIT 


It  is  difficult  to  quantify  what  is  important  to  con- 
tainment of  the  cavity  gases  when  comparing  the  calculations 
of  Table  5 (Section  4).  Some  possible  "figures  of  merit"  are 
discussed  in  this  section.  In  Reference  4,  and  in  Section  4 
of  this  report,  the  values  of  the  peak  residual  stresses 
(particularly  the  hoop  stresses) , relative  to  the  scalar 
overburden  pressure,  were  compared.  It  is  reasonable  to  use 
the  stresses  with  the  overburden  subtracted  out  as  a measure 
of  containment  since  the  philosophy  is  that  it  is  the  addi- 
tional stresses  over  the  overburden  which  forms  the  contain- 
ment membrane  to  prevent  hydrofracture  from  the  cavity.  There 
it  was  stated  that  containment  was  better  when  the  peak 
residual  stresses  were  higher. 

This  simple  figure  of  merit  does  not  give  the  complete 
picture,  however.  Figures  2-7  show  that  the  residual  stress 
fields  look  quite  different  for  strong  media  (where  the  peak 
stresses  are  high)  and  weak  media.  For  weak  media,  the  com- 
pressive hoop  stresses,  though  smaller  in  magnitude,  tend  to 
occur  over  a larger  radial  distance,  i.e.,  the  residual 
stress  membrane  is  thicker.  However,  we  believe  that  a thin- 
ner, but  significantly  higher  stress  membrane  is  far  more  ef- 
fective as  a containment  membrane,  so  that  we  will  continue 
to  use  the  peak  residual  stress  in  our  figure  of  merit. 

Of  course,  these  residual  stresses  must  resist  a cavity 
pressure.  Based  on  the  simple  mechanics  of  materials  problem 
of  a thin  walled  pressure  vessel,  the  ratio  of  cavity  pressure 
to  peak  transverse  residual  stress  would  appear  to  be  a 
reasonable  figure  of  merit.  However,  the  cavity  pressure  com- 
puted using  SKIPPER  is  not  a good  measure  of  the  real  cavity 
pressure.  It  takes  no  account  of  water  getting  into  the 
cavity  from  the  surrounding  rock  or  of  additional  thermal 
equilibration  expected  in  the  rock  outside  the  calculated 
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cavity.  A procedure  developed  by  Peterson  and  Lie^^  to  ob- 
tain cavity  pressures  and  temperatures  for  cavity  leak  cal- 
culations was  modified  here  and  u^ed  to  calculate  a better 
value  of  cavity  pressure  from  the  output  of  the  SKIPPER  cal- 
culations. 

Starting  from  the  final  state  obtained  from  the 
SKIPPER  output  (cavity  radius,  cavity  energy,  internal  energy 
and  mass  of  adjacent  zones,  etc.)  and  beginning  at  the  cavity 
boundary,  small  amounts  of  mass  of  rock  (and  its  accompanying 
internal  energy)  are  equilibrated  with  the  melted  cavity 
material  which  is  at  a specific  interna]  energy  greater  than 
the  melt  energy  of  the  rock.  This  calculational  procedure 
approximates  the  physical  processes  of  spall  and  equilibration. 
The  process  is  stopped  when  sufficient  mass  has  been  added 
so  that  the  equilibrated  "cavity"  is  at  the  melt  energy. 
Finally,  the  volume  of  steam  in  the  new  cavity  is  computed 
based  on  the  water  content  of  the  undisturbed  rock,  and  a 
final  cavity  pressure  is  calculated  (for  the  steam) . 

One  question  which  arises  involves  how  much  water  is 
in  the  undisturbed  rock.  The  water  content  by  weight,  f , 
given  in  Table  2 of  Section  2,  is  defined  as  the  amount  of 
water  that  can  be  removed  at  105°C.  This  does  not  include 
the  bound  water  in  the  rock.  For  Rainier  Mesa  tuff,  this 
amounts  to  approximately  7 percent  extra  water  by  weight. 

The  cavity  pressure  was  computed  for  each  of  the  calculations 
shown  in  Table  5,  using  f alone  and  also  including  the  ad- 
ditional bound  water  (assumed  to  be  7 percent  for  all  media) . 
Table  6 presents  the  results  of  these  computations.  Note  that 
the  effect  of  equilibrating  the  cavity  with  the  outside  mate- 
rial is  to  reduce  the  cavity  pressure  for  all  cases.  When  the 
bound  water  is  included  in  the  computation  of  cavity  pressure, 
less  of  a change  results,  as  would  be  expected.  (The  mass  of 
rock  material  reduces  the  pressure  while  the  water  tends  to 
increase  it.) 


(weak)  0.22  77  44.49  38.5  24.0  32.4  0.50  0.31  0.42 


Here,  ofc  is  the  peak  transverse  residual  stress  and 
the  figures  of  merit  chosen  are  ratios  of  the  calculated 
cavity  pressures,  P^,  and  P^  to  a^.  Let  us  define  a higher 
ratio  of  cavity  pressure  to  as  a greater  relative  contain- 
ment risk.  For  both  the  figures  of  merit  in  which  cavity 
pressure  from  an  equilibrated  cavity  was  used,  both  Baneberry 
calculations  are  clearly  less  conservative  and  therefore 
relatively  greater  containment  risks.  Using  the  unequilibrated 
cavity  pressure,  the  "weak",  best  guess,  Baneberry  2 calcula- 
tion is  the  worst  risk.  However,  the  "strong"  Baneberry  1 
calculation  lies  within  the  data  for  the  other  events. 

Summarizing,  for  all  figures  of  merit  investigated, 
the  "best  guess"  Baneberry  2 calculation  is  clearly  the  worst 
risk.  Only  when  the  unequilibrated  cavity  pressure  (a  hardly 
justifiable  choice)  is  used,  does  even  the  strong  upper  limit 
Baneberry  1 calculation  appear  as  safe  as  for  other  events. 
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